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Abstract 

With the quantum diffusion approach the unexpected behavior of fusion cross section, angular 
momentum, and astrophysical S-factor at sub-barrier energies has been revealed. Out of the region 
of short-range nuclear interaction and action of friction at turning point the decrease rate of the 
cross section under the barrier becomes smaller. The calculated results for the reactions with 
spherical nuclei are in a good agreement with the existing experimental data. 
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I. INTRODUCTION 



The fusion cross section at low energies crucially depends on the capture probability 
of the projectile by the target nucleus, i.e. on the probability to pass the Coulomb bar- 
rier. There are many experimental and theoretical studies of heavy ion fusion reactions 
at extreme sub-barrier energies. The data obtained are of interest for solving the astro- 
physical problems related to nuclear synthesis. Resent measurements of the fusion cross 
sections a at energies E c , m , below the Coulomb barrier jl-4| showed the very rapid fall of 
a just below the barrier. In the ^-factor representation js], S = E cm aexp(2nr]) where 
r/(E c m ) = ZiZ 2 e 2 / {2h? E c m ) is the Sommerfeld parameter, the steep fall-off of the cross 
sections is related to a maximum of the S'-factor. The indications for this maximum have 
been found in Refs. However, its origin is still discussed. The so-called logarithmic 
derivative, L(E cia ) = d(\n(aE cra )) / dE cm , shows a growth at -E c . m . corresponding to the 
maximum of S-factor. 

The experiments jfiS] with the reactions 16 O, 22 Ne+ 208 Pb, where the fusion and capture 
cross sections coincide, demonstrated the decreasing rate of fall of the cross sections at 
energies about 3-4 MeV below the Coulomb barrier. Although this finding has to be checked 
in other experiments, it can not be presently ignored and deserves the theoretical analysis. 
In this paper we will show that unexpected behavior of sub-barrier fusion can be related to 
the switching off the nuclear interaction at external turning point r ex . If the colliding nuclei 
approach the distance Ri n t between their centers, the nuclear forces start to act in addition 
to the Coulomb interaction. Thus, at R < Ri nt the relative motion may be more coupled 
with other degrees of freedom. At R > Ri n t the relative motion is almost independent of 
the internal degrees of freedom. Depending on whether the value of r ex is larger or smaller 
than interaction radius Rint-, the impact of coupling with other degrees of freedom upon the 
barrier passage seems to be different. 

To clarify the behavior of capture cross sections at sub-barrier energies, the further devel- 
opment of the theoretical methods is required. The conventional coupled-channel approach 
with realistic set of parameters is not able to describe the capture cross sections either below 
or above the Coulomb barrier {3]] . The use of quite shallow nucleus- nucleus potential [9| with 
adjusted repulsive core considerably improves the agreement between the calculated and ex- 
perimental data. Besides the coupling with collective excitations, the dissipation, which is 
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simulated by an imaginary potential in Ref. jsj or by damping in each channel in Ref. 10], 
seems to be important. 

The quantum diffusion approach ll| is based on the quantum master-equation for the 
reduced density matrix and takes into account the fluctuation and dissipation effects in 
collisions of heavy ions which model the coupling with various channels. As demonstrated 
in Ref. [Ill ], this approach is successful for describing the capture cross sections at ener- 



gies near the Coulomb barrier. Here, we apply it at wide energy interval including the 
extreme sub-barrier region. To avoid the effects of nuclear deformation, we treat reactions 
16 O, 22 Ne, 48 Ca+ 208 Pb. The collisions of these spherical nuclei is treated in terms of a single 
collective variable: the relative distance R between the colliding nuclei and the conjugate 
momentum P. The coupling with other degrees of freedom is taken into account in the 
friction and diffusion. 



II. MODEL 



We use the nucleus-nucleus potential 



12] which naturally contains the repulsive part 



because of a density-dependent nucleon-nucleon force. The nucleon densities of the projectile 
and target nuclei are specified in the form of the Woods-Saxon parameterization, where the 
nuclear radius parameter is tq = 1.15 fm and the diffuseness parameter takes the values 
a = 0.55 fm for the 208 Pb, 48 Ca, 22 Ne nuclei and a = 0.53 fm for the 16 nucleus. The 
nucleus-nucleus potential V has quite a shallow pocket (Fig. 1). As the centrifugal part 
of the potential grows, the pocket depth becomes smaller, while the position of the pocket 
minimum moves towards the barrier at R = Rh. This pocket is washed out at angular 
momenta J > 80. As demonstrated in Refs.p, [3], the internuclear potentials show a 
significant deviation from the conventional Woods-Saxon shape. The thicker the potential 
is, the smaller the penetrability is, and also the stronger the energy dependence of the 
penetrability is. The thick potential barrier obtained for the 16 O+ 208 Pb reaction is thus 
consistent with the recent experimental observations [9] that the fusion excitation function 
is much steeper than theoretical predictions at deep sub-barrier energies. The present study 
suggests that the origin of the steep fall-off of fusion cross section is partly attributed to the 
deviation of the internuclear potential from the Woods-Saxon shape. 
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The capture cross section is a sum of partial capture cross sections 



&c{E c .m.) = 0"c (-5c.nl., J) 



7rA 2 ^(2J+l)P cap (P c . m ., J), 
j 



(1) 



where A 2 = h 2 / (2jiE c-m ) is the reduced de Broglie wavelength and the summation is in 
possible values of angular momentum J at given bombarding energy £? c .m.- The partial 
capture probability P cap is defined by the passing probability of the potential barrier in the 
relative distance P at given J. 

The value of P cap can be obtained by integrating the propagator G from the initial state 
(Rq, Pq) at time t = to the final state (R, P) at time t: 



P 



cap 



lim 

t— >oo 



dR / dP G(R,P,t\R o ,P o ,0) 



lim -erfc 

t— >oo 2 



R(t) 



line 



in 



is 



obtained by using the propagator G 



(2) 



The second 

7r- 1 |detS- 1 | 1/2ex P(-q Ts_1 q) (QR.(t) = R-R^j, q P (t) = P-P({), R(t = 0) = Rq, 
P(t = 0) = P , £y(t) = 2q i (t)q j (t), E y (t = 0) = 0, i, j = P, P) calculated in Ref. U for 
the inverted oscillator which approximates the nucleus-nucleus potential V in the variable 
P. The frequency u> of this oscillator with internal turning point r m is defined from the 
condition of equality of the classical actions of approximated and realistic potential barriers 
of the same hight at given J. It should be noted that the passage through the Coulomb 
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4l7j. This 



barrier approximated by a parabola has been previously studied in Refs. 
approximation seems to be well justified for the reactions considered. 

All quantum-mechanical, dissipative effects and non-Markovian effects accompanying the 



passage through the potential barrier are taken into consideration in our formalism . The 
generalized fluctuation-dissipation relations contain the influence of quantum effects on the 
collective motion. We address the dynamics of damped non-Markovian collective mode in 



terms of the first moment R(t) (the average of the coordinate P) and variance £##(£) in the 
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coordinate: 



R(t) 



A t R + B t P 

t 

drB> 



2fr 2 A 7 2 



B f 



71 

x coth 
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dr B, 



dn 



n 



fi 2 + 7 2 



o 

2T 



cos[Q(t — T )], 



p t=l 

3 

^A[si(s i + 7) + ^A 7 ///]e Si *. 



(3) 



.0 



The derivation of equations for R(t) and is presented in Refs. 1 161 ] as well as 

in Appendix A. Here, R(0) = R , Err{0) = 0, A = 1, and £ = 0. In Eqs.©, ft = 
[(si - s 2 )(si - S3)]" 1 , ft = [(s2 - si)(s 2 - s 3 )] _1 and /3 3 = [(s 3 - si)(s 3 - s 2 )] _1 , and s; are 
the real roots (si > > s 2 > S3) of the following equation: 

(5 + 7) (s 2 - to 2 ) + h\-fs/fx = 0. (4) 

Here, /1 is the reduced mass, u 2 . = u 2 {l — hXy/\jj,(si + 7 )(s 2 + 7)]} is the renormalized 
frequency in the Markovian limit, and the value of A is related to the strength of linear 
coupling in coordinates between collective and internal subsystems. The friction coefficient 
in R is set as HX = —h(s\ + s 2 ) = 2 MeV. It has the value close to those calculated within 
other approaches [18j. Because of the quite large values of u, the calculated results are 
rather insensitive to the value of temperature T < 1 MeV. Non-Markovian effects appear in 
the calculations through the internal-excitation width hy = 15 MeV. Since the relaxation 
time for the internal subsystem is much shorter than the characteristic time of collective 
motion, the condition 7 ^> 00 should be fulfilled. The value of A is partly related to the 
value of 7. In the limit of A — > the value of 7 should go to infinity to meet the Markovian 
dynamics. 



The Eqs. © and © result 
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(5) 



where i/j(z) is digamma function. Since Err — > oo at t — > oo, the value of P cap is independent 
of r; n . Using Eq. fll]), in the limit of small temperature (T — Y 0), which is suitable for sub- 
barrier fusion, we obtain 



P> cap 2 *^f < -' 



TTSl (7 - Si) \ V2 ^0 R o/ S l + P 



2/^K 2 -^i)/ [ 7 ln( 7 M)] 1/2 



(6) 



This expression is used to calculate P cap in our paper. The conventional approach to mod- 
eling fusion cross sections is based on calculating quantum mechanical barrier transmission 
probabilities for each partial wave of relative motion, without regard (in first approxima- 
tion) to coupling to internal degrees of freedom. If the coupling with internal degrees of 
freedom is disregarded in the quantum diffusion approach at zero temperature, then A — > 0, 
si — > uq, and 7 — > 00 so that 4hXln( , ~f /uq) — > 1. In this case the well-known expression 
-Pcap ~ exp[— 27r(Vf> — E c . m ) / (Jtojq)] is obtained where V& is the hight of the barrier at given J. 
Thus, our diffusion approach contains the quantum mechanical barrier transmission proba- 
bility. In the presence of coupling with internal degrees of freedom the capture at sub-barrier 
energies also occurs due to the quantum noise. The actions of friction and diffusion are the 
reverse of each other ll| in the dissipation of energy. As a result, in the reactions treated 



the dissipation of energy is negligible at — £ c . m . > 3.5 MeV. 

The nuclear forces start to play a role at Ri n t = Rb + 1-1 fm where the nucleon density of 
colliding nuclei approximately reaches 10% of saturation density. In Fig. 2 the interaction 
radius Ri nt as well as internal and external turning points are shown for the 16 O+ 208 Pb 
reaction at zero angular momentum. If the value of r ex corresponding to external turning 
point is larger than interaction radius Ri n t, we take Ro = r ex and Pq = in Eq. <Q. For 
r ex < Rint, it is naturally to start our treatment with Rq = Ri nt and Pq defined by the kinetic 
energy at R = Rq. In this case the friction hinders the classical motion towards smaller R. 
If Pq — at Rq > Rint, the friction almost does not play a role in the transition through 
the barrier. Thus, two regimes of interaction at sub-barrier energies differ by the action of 
nuclear forces and the role of friction at R = r Px . 



III. CALCULATED RESULTS 



Besides the parameters related to the nucleus- nucleus potential, two parameters Try = 15 
MeV and h\ = 2 MeV are used for calculating the capture probability. All calculated results 
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are obtained with the same set of parameters and are rather insensitive to the reasonable 
variation of them [ill , [if]]. 

In Figs. 3-5 the calculated capture cross sections for the reactions 16 22 Ne, 48 Ca+ 208 Pb 



19r|2l|. There 



are in a good agreement with the available experimental data 0,0, 6-8, 
is sharp fall-off of the cross sections just under the barrier. With decreasing E c _ m _ up to 
about 3.5-5.0 MeV below the Coulomb barrier the regime of interaction is changed because 
at external turning point the colliding nuclei do not reach the region of nuclear interaction 
where the friction plays a role. As a result, at smaller E cm the cross sections fall with 
smaller rate. With larger value of Ri n t the change of fall rate would occur at smaller E c , m _. 
However, the uncertainty in definition of Ri nt is rather small. Therefore, the effect of the 
change of fall rate of sub-barrier fusion cross section should be in the data if we believe that 
friction start to act only when the colliding nuclei approach the barrier. Note that at energies 
of 5 MeV below the barrier the experimental data have still large uncertainties to make a 
firm experimental conclusion about this effect. The effect seems to be more pronounced in 
the collisions of spherical nuclei. The collisions of deformed nuclei occurs at various mutual 
orientations from which the value of Ri nt depends. 

The calculated average angular momenta (J 2 ) = Yl j J{J + ^) a c{.E cia i J) / a c {E cia ) of 
fused systems versus E c _ m , are presented in Figs. 3-5 as well. At energies of 3-4.5 MeV 
below the barrier (J 2 ) has a minimum. The experimental data indicate the presence 
of the minimum as well. On the left-hand side of this minimum the dependence of (J 2 ) of 
i?c.m. is rather weak. The similar weak dependence has been found in Ref. 23J at extreme 
sub-barrier region. Note that the found behavior of (J 2 ), which is related to the change of 
the regime of interaction between the colliding nuclei, would affect the angular anisotropy 
of the products of fission following fusion. 

In Fig. 6 the functions L(E cra ) and S(E cra ), and the fusion barrier distribution 
d 2 (E c _ m ,a) j dE\ m are presented for the 16 O+ 208 Pb reaction. The logarithmic derivative 
strongly increases just below the barrier and then has a maximum. This leads to the maxi- 
mum of ^-factor which is seen in the experiments j^]. After this maximum ^-factor decreases 



with E C Ta . and then starts to increase. The same behavior has been revealed in Refs. [24] 
by extracting the S'-factor from the experimental data. If for finding the logarithmic deriva- 
tive we use only the cross sections calculated at the energies where the experimental cross 
sections are available, the function L(E c _ m _) would be similar to that obtained with the ex- 
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perimental data [3]. Therefore, the energy increment in Ref. {^J seems to be too large to 
extract a function L(E c . m .) with a very narrow maximum. This increment should be at least 
0.2 MeV. The fusion barrier distribution calculated with this small energy increment has 
only one maximum. Using larger energy increment of 0.6 MeV, one can get few oscillations 
in <P(E c . m a)/dEl m , 



IV. SUMMARY 



The quantum diffusion approach has been applied to study the capture or fusion cross 
sections at sub-barrier energies. The available experimental data at energies above and below 
the Coulomb barrier are well described. Due to the change of the regime of interaction (the 
turning-off the nuclear forces and friction) at sub-barrier energies, the decrease rate of the 
cross sections is changed at about 3.5-5.0 MeV below the barrier. This change is reflected 
in the functions (J 2 ), L(E c . m ), and S(E c _ m _). The average angular momentum of compound 
nucleus versus -E c .m. would have a minimum and then saturation at sub-barrier energies. 
This behavior of (J 2 ) would increase the expected anisotropy of angular distribution of the 
products of fission following fusion. The energy increment of 0.2 MeV has to be used in the 
experiment to get the cross sections suitable for calculating the value of L and the barrier 
distribution. 

This work was supported by DFG and RFBR. The Polish- JINR and IN2P3-JINR Coop- 
eration programs are gratefully acknowledged. 



Appendix A 



For quantum nuclear system, the Hamiltonian H depending explicitly on the collective 



coordinate R, the canonically conjugate co 
of freedom has been constructed in Refs. 



lective momentum P, and on the internal degrees 
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161 ] . Using this Hamiltonian, one can obtain 



the system of integro-differential stochastic equations for the Heisenberg operators R and P 

P(t) 



R(t) 



P(t) = (iu%R(t) - / drK(t - t)R(t) + F(t). 



(Al) 
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Here, K(t) = h\'ye~' yt is a dissipation kernel and F(t) is a random force. This system of 
Eqs. (Al) is the system of generalized nonlinear Langevin equations. The integral term 
in the equations of motion means that the system is non-Markovian and has a "memory" 
of the motion in the trajectory preceding the instant t. The operator of the random force 
F has the form of a Gaussian distribution with a zero mean, ((F(t))) = 0, and a nonzero 
variance, 



~ oo 

f ^ 



«F(t)F(0» = / dQ——^ .coth 



7T 







fi 2 + 7 2 



2T 



cos[ft(t -*')]. (A2) 



The symbol ((•)) denotes the mean over the internal degrees of freedom. To solve Eqs. 
(Al) analytically, we use the Laplace transform method. After finding expressions for the 
images, we obtain explicit expressions for the originals Ll[ Q , 

t 

R(t) = A t R(0) + B t P(0) + J drB T F(t - r), (A3) 

Q 

where the explicit expressions for A t and B t are given in Eq.(3). Using the time dependence 
R(t), we obtain the values R(t) (see Eq.(3)) and 

t t 

E RR (t) = 2(R(t) - R^)y = 2 J dr'B T> j dr" B T „ ((F(t - r')F(t - /))) (A4) 

o o 

averaged over the whole system. Using Eq.(A2), Eq.(A4) can be rewritten in the form of 
Eq.(3). 
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FIG. 1: The nucleus- nucleus potentials calculated at J = (solid curve), 30 (dashed curve), 60 
(dotted curve), and 90 (dash-dotted curve) for the 16 O+ 208 Pb reaction. 

FIG. 2: The nucleus-nucleus potential calculated at J = (solid lines) for the 16 O+ 208 Pb reaction. 
The position Ri, of the Coulomb barrier, radius of interaction Ri n t, and external and internal 
turning points for some value of E c _ m , are indicated. 

FIG. 4: The same as in Fig. 3, but for the 22 Ne+ 208 Pb reaction. The experimental data marked 
by closed squares and circles are from the two runs of Ref . [8fl . 



FIG. 5: The same as in Fig. 3, but for the 48 Ca+ 208 Pb reaction. T 
by closed circles and squares, and open circles are taken from Refs. 



re experimental data marked 



19H2l|, respectively. 



FIG. 6: The calculated values (solid lines) of logarithmic derivative L (upper part), astrophysical 
S'-factor (middle part) with t]q = r](E c . m , = Vf,), and fusion barrier distribution (lower part) for 
the 16 O+ 208 Pb reaction. The dashed line shows^the values of L obtained only with the calculated 
cross sections at E c m _ used in the experiment 



FIG. 3: The calculated (solid lines) fusion cross section versus S c .m. (upper part) and average 
angular momenta of compound nucleus (lower part) versus -B c . m . for the 16 O+ 208 Pb reaction are 
compared with the experimental data. The experimental cross sections marked by open squares and 



circles, and closed rhombus and triangles are from Refs. 3j, |4, 
values of (J 2 ) (solid squares) are from Ref. 
by arrow. 



I, 111 , respectively. The experimental 



221 ] . The value of the Coulomb barrier Vb is indicated 
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